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Scope and Method of Study: The performance of an axial vane- type 
swirler was investigated to aid in computer modeling of gas tur- 
bine combustor flowfields and in evaluation of turbulence models 
for swirling confined jet flow. The swirler studied is annular 
with a hub-to-swirler diameter ratio of 0,25 and ten adjustable 
vanes of pitch-to-chord ratio 0,68. Measurements of time-mean 
axial, radial, and tangential velocities were made at the swirler 
exit plane using a five-hole pitot probe technique with computer 
data reduction. Nondi mens ionali zed velocities from both radial 
and azimuthal traverses are tabulated and plotted for a range of 
swirl vane angles <}) from 0 to 70 degrees. In addition, a study 
was done of idealized exit-plane velocity profiles relating the 
swirl numbers S and s' to the ratio of maximum swirl and axial 
velocities for each idealized case, and comparing the idealized 
swirl numbers with ones calculated from measured profiles. 

Findings and Conclusions: Measurements of time-mean velocity components 
at the swirler exit plane show clearly the effects of centrifugal 
forces, recirculation zones, and blade wakes on the exit-plane 
velocity profiles. Assumptions of flat axial and swirl profiles 
are found to be progressively less realistic as the swirl vane 
angle increases, with axial and swirl velocities peaking strongly 
at the outer edges of the swirler exit and significant non-zero 
radial velocities present. Higher-order idealized profiles gave 
improved correspondence with moderate tp high swirl cases, but none 
of the idealizations studied could approximate the measured pro- 
files satisfactorily. For strong swirl, the central recirculation 
zone extended upstream of the exit plane, and nonaxi synmetry was 
found in all swirl cases investigated. 
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MOMENCLATURE 
English Symbols 

blade chord width 
swirl er exit diameter 
test section diameter 
velocity ratio w^/Uq for case I 

axial flux of momentum; velocity ratio w^^/u^ for case II 

W-/u for cases III V 
mov mo ^ ■ ■" 

time-mean pressure, N/ni = Pa 
blade spacing or pitch 
swirl number = Gg/(G^d/2) 

axial, radial and tangential components of velocity 
axial, radial, azimuthal cylindrical polar coordinates 
hub-to-swirler diameter ratio d^^/d 

Greek Symbols 

yaw angle of probe - tan“^ (w/u) 
pitch angle of probe = tan"^ Cv/(u^ + 
azimuth angle 
density 

pitch - to - chord ratio 

swirl vane angle - tan (w^^/U|^), assuming perfect vanes 
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Subscripts 

atm ambient atmospheric conditions 

C,N,S,E,W center, north, south, east, west pitot pressure ports 

h hub 

in inlet conditions, upstream of swirl er 

m maximum profile value 

0 value at swirl er outlet 

X axial direction 

0 tangential direction 

« reference value at edge of swirl er exit 

Superscripts 

' alternate form, neglecting pressure variation; 

fluctuating quantity 

— time-mean quantity 
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CHAPTER I 
INTRODUCTION 

1.1 Combustor Flowfield Investigations 

The problem of optimizing gas turbine combustion .chamber design is 
complex, because of the many conflicting design requirements. The need 
for a more complete understanding of the fluid dynamics of the flow in 
such combustion chambers has been recognized by designers in recent 
years, and research is continuing on several fronts to alleviate the 
problem. 

As part of an on-going project at Oklahoma State University, studies 
are in progress concerned with experimental and theoretical research in 
2-D axi symmetric geometries under low speed, nonreacting, turbulent, 
swirling flow conditions. The flow enters the test section and proceeds 
into a larger chamber (expansion ratio D/d = 2) via a sudden or grad- 
ual expansion (side-wall angle a = 90 and 45 degrees). Inlet swirl vanes 
are adjustable to a variety of vane angles with = 0, 38, 45, 60 and 
70 degrees being emphasized. The general aim of the entire study is to 
characterize the time-mean and turbulence flowfield, recommend appro- 
priate turbulence model advances, and implement and exhibit results of 
flowfield predictions. The present contribution concentrates on the 
time-mean flow characteristics being generated by the upstream annular 
swirler, using a five-hole pitot probe technique. 
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1.2 Previous Studies 


ORIGIN/'iL PmE m 
OF POOR QUAUTY 


Research is progressing in several areas related to the flow facil- 
ity investigation just described. Computer simulation techniques are 
being used to study the effect of geometry and other parameter changes 
on the flowfield. An advanced computer code (1) has been developed to 
predict confined swirling flows corresponding to those studied experi- 
mentally. Tentative predictions (2) have now been supplemented by pre- 
dictions made from realistic inlet conditions (3) for a complete range 
of Swirl strengths with downstream nozzle effects (4), Accuracy of pre- 
dictions from a computer model is strongly dependent on the inlet boun- 
dary conditions used, which are primarily determined by the swirl er and 

its performance at different vane angle settings. In the earlier predic 
tions, the velocity boundary conditions at the inlet to the model com- 
bustor were approximated by idealized flat profiles for axial and swirl 

velocity, with radial velocity assumed to be zero. However, recent 
measurements taken closer to the swirler exit show that the profiles 
produced are quite nonuniform, with nonzero radial velocity and nonaxi- 
symmetry . 

The flowfield in the test section is being characterized experimen- 
tally in a variety of ways. Flow visualization has been achieved via 
still (5) and movie (6) photography of neutrally buoyant helium-filled 
soap bubbles and smoke produced by an injector and a smoke wire. Time- 
mean velocities have been measured with a five-hole pitot probe at low 
(5) and high (7) swirl strengths. To help in turbulence modeling, com- 
plete turbulence measurements have been made on weakly (8) and strongly 
(9) swirling flows, using a six-orientation single-wire hot-wire tech- 
nique. An alternative three-wire technique has also been shown to be 
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useful in the eoniplex flow situations (IQ). 

References to previous work done elsewhere are found in Chapter II, 
relating to theoretical analysis of swirler performance. 

1*3 Scope and Objectives 

A key element in swirling flow studies is the swirl generator used. 

Since it lies at the inlet to the combustor model, the swirler can have 
a strong influence on the measurements or predictions made downstream, 

Better definition of the swirler's performance characteristis^s is needed. 

In the present study, the main objective has been to make time-mean 
velocity measurements as close as possible to the swirler exit, so as to 
define more accurately the performance characteristics of the swirler, 

A range of swirl-blade angles f|t from 0 to 70 deg, is considered. Speci- 
fic objectives include: 

1. Investigate the flow turning effectiveness of flat blades in 
annular axial vane swirl era at various blade angles, cj), 

2. Investigate the degree of iionaxi symmetry introduced by vane- 

type swirl ers. | 

3. Establish correlations between the blade angle «f) and the velo- 
city profiles and degree of swirl actually produced. 

4. Evaluate the applicability of idealized velocity profiles used 
recently in flowfield prediction codes, and specify more rea- 
listic idealized profiles for future use. 

5. Provide swirler exit data usable as inlet conditions in predic- 
tion codes being used to establish, evaluate, and improve 
turbulence models* 
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1,4 Out1ir> of the Thesis 

In the previous sections, the scope and objectives of this study 
were presented, with the significance of the study in relation' to past 
and present work on combustor flowfield investigations being high- 
lighted. 

Chapter II describes nathematical derivations from Idealized 
swirler exit velocity profiles, relating the swirl number to the ratio 
of maximum swirl and axial velocities for several cases. 

Chapter III covers the experimental equipment and procedures used 
for measurement of the swirler exit flowfield. It includes descrip- 
tions of the flowfield facility, the swirler, and the five-hole pitot 
probe and its associated instrumentation. Calibration, measurement, 
and data reduction procedures are also briefly described. 

The first two sections of Chapter IV discuss experimental results 
from radial and azimuthal traverses, respectively, noting the presence 
of nonaxi symmetry, recirculation, and strong velocity gradients at the 
swirler exit plane. A third section describes the results of a check 
on sensitivity of the measurements to calibration errors. The last 
section of Chapter IV compares the swirl numbers calculated from mea- 
sured profiles and from the idealizations of Chapter II to judge the 
usefulness of the idealized profiles. 

Chapter V presents conclusions drawn from the above results and 
makes recommendations for further research on this topic. 

Appendixes A and B include tables and figures, respectively. A 
description of revisions to be computer program for reduction of five- 
hole pitot probe data is in Appendix C, and a listing of the program 
with sample input is in Appendix D. 
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CHAPTER II 

IDEALIZED PROFILE DERIVATIONS 
2,1 Idealized Velocity Profiles 

All theoretical analyses of sv/irler performance and most numerical 
simulations of combustor flowfields have used simple idealized swirler 
exit velocity profiles. Common assumptions made include flat axial and 
swirl velocity profiles downstream of the swirler for swirlers with vanes 
of constant angle (2, 5, 11, IH), and flat axial profile with linear 
swirl profile (solid-body rotation) for swirlers with helicoidal vanes 
and for tangential -entry swirl generators (13, 14) , These, however, have 
been shown to be qi'ite unrealistic (3, 12, 15) and to lead to consider- 
able errors in computer siitulations (4). Although the best approach for 
numerical simulations is to use experimentally measured profiles if they 
are available, idealized profiles are very useful in theoretical work. 

If more realistic profile t.ssumptions can be developed which are still 
mathematically tractable, niore useful analytical results may be derived. 
Better idealized profiles would also be useful as inlet boundary condi- 
tions for computer modeling when measured data is not available. 

Measurements have shown (3) that linear and parabolic profiles of 
axial velocity are more appropriate for moderate and high swirl cases, 
and that the swirl velocity also approaches a parabolic profile at high 
swirl Strengths, with most of the flow leaving near the outer boundary 
of the swirler. Several combinations of linear and parabolic idealized 
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profiles are shown in Figure 1, along with the flat and linear profile 
assumptions used in previous studies- Parameters associated with these 
profiles are investigated in Section 2,3. 

2,2 Definition of Swirl Parameters 

The swirl number is a nondimenslpQal parameter used to characterize 

the degree of swirl generated by a swirl er. It is defined as follows 

(13): G. 

S = — 

( 1 ) 

where the axial flux of angular momentum Gg is given by 

27t d/2 . 

Gg - / de / [^luw + puV] r ' dr (2) 

and the axial flux of axial momentum G^ is given by 

2tt d/2 « — ip 

= / dD * [,)U + + (p - P^) ]r dr (3) 

and d/2 is the swirler exit radius (4). These equations are obtained 
from appropriate manipulation of the axial and azimuthal momentum equa- 
tions, respectively. In free jet flows these two expressions are 
invariant with respect to downstream location. In the axial momentum 
expression, the pressure term (p - p^) is given from radial integration 
of the radial momentum equation (16) by 

(P - Pj = /'" [pw^ 1] dr - pv'^ ... 

d/2 

If the pressure term is omitted from the axial momentum, the dynamic 
axial momentum flux 6 is obtained: 

A 

I 2tT q/ 2 p n 

Gy = / de / + pu'‘^] r dr 

^ 0 


( 5 ) 
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This leads to an alternate definition of swirl number (17): 
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" G;(d/2)' (6) 

If turbulent stress terms are neglected, it is apparent that a know- 
ledge of the distribution of the time-mean u and w velocity components 
across the swirler is sufficient to calculate either swirl number. The 
idealized exit velocity profiles provide just such knowledge, and expres- 
sions relating swirl number to the ratio of maximum exit swirl and axial 
velocities can now be derived for each of the profile types. As the 
procedure is similar for each of the five cases, a detailed derivation 
will be shown for the first case only, with only final results given 
for the other four. 

2.3 Swirl Numbers for Idealized Profiles 

By assuming axisymmetric flow and neglecting turbulent stresses 
as stated previously, the definitions in Equations (2) through (4) re- 
duce to 


c— 

CM 

II 

O 

CD 

d/2 2 

[ouw]r dr 

') 

(7) 


d/2 2 

/ [nu •!* (p - p^)]r dr 

(8) 

(P - Pj 

- f [i>w^ i]dr 

(9) 


d/2 


When the expressions for axial and swirl velocity for case I (see Figure 
1) are substituted into Equation (7), one obtains 
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Gq = fitpu^w (d/2)'^ (10) 


Substitution of w(r) - w. into Equation (9) and integrating produces 

0 

(p - pj = pwf [ln(f) - ln(d/2)] (11) 


After substituting Equation (11) into Equation (8) and integrating, the 
expression becomes 

= itpu^ (d/2)2 [1 - |- (^) ] (12) 

0 

Finally, putting Equations (10) and (12) into Equation (1) and defining 
the velocity ratio F = w^/u^, the swirl number S can be expressed thus: 

. _ 2F/3 

1 - rn 

The alternate swirl number S' follows from finding the dynamic axial 
flux of axial momentum: 

Gj^ = ITPU^ (d/2)^ (14) 

Using this in Equation (6) leads to the simple expression, 

S’ = 2F/3 (15) 


By the same procedure, expressions for S and S' for the other four 
cases are found to be as follows: 


For case II with u(r) = u^, w(r) = defining G as 


w /u : 
mo' 0 


S - 


G/2 


(16) 


1 - G74 
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and 


S* s G/2 
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07) 


For case III with u(r) = i^^). w(r) = 

, and defining 

” V^V= 


c _ 4H/5 

1 - H^/2 

(18) 

and 


S'= 4H/5 

(19) 

r 2 

For case IV with u(r) = w(r) - . 

and defining I 

Vmo = 


c . I 

b - -y 

1 - 3r/4 

(20) 

and 


S' = I 

(21) 

r 2 

Finally, for case V with u(r) = u^^^ {"^^) , w(r) = 

r 2 
”mo ’ 

and defining J as 


. _ 40/7 

5 P 

1 - 20 V3 

(22) 

and 


S' = 40/7 

(23) 


Each of these expressions for S and S' may be inverted to yield 
the velocity ratio as a function of swirl number. A summary of the 
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Case I - 

c "4/ (3S) [ 4 /(3S) l^ + 8 

2 

F « 3S'/2 

Case II - 

G = -2/(S) + [Z/(S)f- t 16 

t, 2 

G = 2S' 


Case III - 

u - -8/(5S) [8/(5S)3^ ^ 8 

H - 2 


H = 5S'/4 


Case IV - 

T - -4/(3S) j- [4/(3$)]^^16/3 
^ " 2 

I = S' 

Case V - 

1 - -6/ (7S) -f C6/(7S)]^ + 6 
J - 2 


J = 7S74 


n 


Numerical values from each of these expressions are given In Table I, 
and the same relationships are shown graphically In Fig. 2 for a range 
of commonly-encountered swirl numbers. 

It Is evident from the equations alone that the S' expressions are 
all simple linear relations. The parameters F through J will Increase 
without bound as the swirl number Is Increased in each case. In con- 
trast, the parameter variation with S shows asymptotic behavior; the 
exit velocity ratios all approach definite values as swirl number In- 
creases. The asymptotic values are also given in Table I. 

Although the curves are generally similar in shape, some observa- 
tions can be made. The curves for cases II and IV are the upper and 
lower extremes for both the S and S' relations, with the curves for 
cases I, III, and V falling In between. This may be anticipated since 
the w profile Is of higher order than the u profile for case II (that 
is, linear versus constant) and the opposite Is true for case IV 
(linear versus parabolic). In the other three cases the u and w pro- 
files are of the same order. 

In appraising the usefulness of ne Idealized profiles, comparison 
may be made with the measured profiles given later in Chapter IV. As 
the swirl strength increases from 0 to 70 deg., corresponding profiles 
of cases I to V appear roughly appropriate. The moderate swirl case 
((})- 45 deg.) gives the best match With Its corresponding Idealization 
(case III, linear axial and swirl profiles), by visual inspection alone. 
However, the presence of the hub and central recirculation zone prevent 
adequate representation by the Idealized profiles, as demonstrated by 
the experimental results discussed in Chapter IV. 
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EXPERIMENTAL EQUIPMENT AND PROCEDURE 
3.1 Combustor Flowfield Facility 

The installation on which all tests were performed is a low-speed 
wind tunnel designed and built at Oklahoma State University. It pro- 
duces uniform flow of relatively low turbulence intensity, with contin- 
uously adjustable flow rats. The facility consists of a filtered intake, 
an axial blower, a stilling chamber, a turbulence management section, 
and a conteured outlet nozzle, A schematic of the facility is shown in 
Fig. 2. 

The intake consists of a rounded entrance containing fixed inlet 
guide vanes, surrounded by a coarse-mesh screen box covered with foam 
rubber panels to filter the incoming ambient air. The blower is a six- 
bladsd propeller-type fan, driven by a 5 h.p. U.S. Varidrive motor which 
can be continuously varied from 1600 to 3100 rpm. 

Air from the blower is expanded into the stilling chamber and passes 
through several fine mesh screens to help remove the turbulence gener- 
ated by the blower. The turbulence level is further reduced by passage 
through the turbulence management section. This section, a round duct 
of 76 cm diameter, contains a perforated aluminum plate (2 mm diameter 
holes) followed by a fine mesh screen, a section of packed straws 12.7 
cm long, and five more fine mesh screens. Most of the turbulence reduc- 
tion occurs in this section, and any traces of fan-induced swirl are 
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effectively removed by the straws. 

To reduce the duct diameter down to the 15 cm outlet diameter, a 
specially contoured nozzle is used, This was designed after the method 
of Morel (18) to minimize boundary layer growth and produce a uniform 
top-hat profile, with no separation or instabilities upstream. The 
nozzle is of molded fiberglass with a steel flange at the outlet for the 
attachment of various test articles. A 1 cm diameter hole a short dis- 
tance upstream of the outlet allows for insertion of a standard pitot- 
static probe to measure the dynamic pressure upstream of the swirler. 

This measurement, with a small correction for difference In flow area, 
is used to calculate the swirler inlet reference velocity, u^.^. 

3,2 Swirler 

The swirler used in this study is annular with hub and housing dia- 
meters of 3.75 and 15.0 cm respectively, giving a hub-to-swirler dia- 
meter ratio z of 0.25. The hub has a streamlined parabolic nose facing 
upstream and a blunt base (corner radius approximately 2 mm) facing down- 
stream. It is supported by four thin rectangular-section struts or 
spider arms from the housing wall. The base of the hub protrudes approx- 
imately 3 mm downstream of the swirler exit plane. Photographs are sche- 
matics of the swirler are shown in Figures 3 through 5. 

The ten vanes or blades are attached to shafts which pass through 
the housing wall and allow individual adjustment of each blade's angle. 
The standard vanes are wedge-shaped for nearly-constant pitch-to-chord 
ratio a of approximately 0.68, which according to two-dimensional cas- 
cade data should give good flow-turning effectiveness. Sets of 
vanes with chord widths of 0.5 and 0.75 of the standard width may be 
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Installed to study the effect of increased pitch-to-chord ratio on turn- 
ing effectiveness, nonaxi symmetry, and radial secondary flow patterns. 
Vane planforms are shown in Figure 6. 

3.3 Five-Hole Pitot Probe and Instrumentation 

Velocity profile measurements were made using a five-hole pitot 
probe (Model DC-1 25-1 2-CD by United Sensor Division of United Electrical 
Controls Co.), one of the few instruments capable of measuring the mag- 
nitude and direction of the local time-mean velocity vector simultan- 
eously. Detailed explanations of five-hole pitot operating techniques 
and basic principles may be found in Reference 5. A schematic of the 
probe tip geometry showing the velocities and angles measured is given 
in Figure 7. 

The probe is mounted in a traversing mechanism (Model Cl 000-1 2 from 
United Sensor) which in turn is im^unted on a 30-cm diameter plexiglass 
tube which fits closely over the swirler exit flange. This tube com- 
prises the test section for combustor flowfield modeling in related 
studies (1-10) and creates confined-jet conditions downstream of the 
swirler. The presence of the '^est section tube has neglftrlble effect 
on the flow patterns observed at the swirler exit plane. 

The traversing mechanism anows the probe to be translated verti- 
cally (on a radial line outward from the test section axis) and rotated 
360 degrees about the probe's yaw axis. In addition to the motion per- 
mitted by the traverse mechanism, the test section tube on which the 
traverse mechanism is mounted may be rotated about its axis with respect 
to the swirler, thereby allowing azimuthal traverses to be performed. 

Tubing from the probe's five pressure taps is routed through selec- 
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tor valves so that pressure differences between any two of the probe's 
five holes may be measured by a differential pressure transducer (Type 
590 Darocel Pressure Sensor by patametrics Inc,, + 10 torr range). The 
resulting pressure difference values are then read directly from a digi- 
tal voltmeter with selectable averaging time -constant (Model 1076 True 
RMS VoltniOter by TSI, Inc.), 

3,4 Cali Oration i Measure!i)ent, and 
Reduction Procedure 

Calibration of the five-hole probe is done using a small free jet 
which has a contoured nozzle similar to that of the flowfield facility. 
The probe tip is placed in the uniform parallel flow of the jet potential 
core and adjusted to zero yaw angle. The probe is then rotated about 
its pitch axis and values of (p^^ - Pg), (Pq - p|^), and (p^ - p^^^^^) pres- 
sure differences are measured at different values of pitch angle 6. 

Velocity measurements with the five-hole probe are made after the 
probe has been carefully aligned with the facility and the pressure 
transducer properly zeroed. At each measurement location, the probe is 
aligned with the local flow direction in the horizontal plane by nulling 
the pressure difference (pjj - p^^). The value of yaw angle 3 is then 
read from the rotary vernier on the traverse mechanism. Finally, values 
of the pressure differences (p^^ - Pg), (p^^ - p^^), and (p^ " Patm- 
measured. 

The raw pressure data are reduced by a computer program to yield 
nondimensionalized values of the u, v, and w velocity components, as 
well as the static pressure at each location. The reduction program 
also performs numerical integration on the radial traverses to obtain 
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values of the axial and angular momentum fluxes, and from these calcu- 
lates the swirl numbers S and S'. Some details of the reduction proce- 
dure are given In Appendix C, the description of changes made to the 
reduction code, while more general descriptions of the original code 
are found in references (19) and (20), A listing of the code with 
sample input anu output Is given in Appendix D. 



CHAPTER IV 


EXPERIMENTAL RESULTS 

« 

Velocity pi^ofiles from both radial and «;izimuthal traverses for each 
of the flowfields investigated are now presented and discussed. 

Table II gives a summary of the operating conditions used during the 
studies. With nonswirling conditions, the low fan speed delivers rela- 
tively high axial velocity and corresponding Reynolds number. At pro- 
gressively higher swirl strength conditions, progressively higher fan 
speeds are used, but even so exit velocities and Reynolds numbers reduce 
because of increasing flow restriction of the swirl er. However, based 
on a limited study elsewhere (4), it is expected that all flowfields 
are in the Reynolds number independent regime. 

The radial traverses consist of ten points from the centerline to 
the swirl er exit radius, spaced 7.6 mm apart. Of these ten, only seven 
stations were actually measured since the hub blocked the inner three 
positions. The azimuthal traverses contain nine points spaced 6 degrees 
apart at a constant radial distance from the centerline. Azimuth angles 
0 were taken from -24 to +24 degrees, with the 9=0 position in line 
with the shaft of one of the swirl vanes. A diagram showing the tra- 
verse patterns on the face of the swirl er is given in Figure 8. 

Unless otherwise stated all traverses are taken immediately after 
the swirler exit downstream face with no expansion blocks present. 
Nominally, this location is x/D = -0.109, where the positon x/D = 0.0 
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IS the expansion station, separated from the swirler in practice (S^-IO) 
with one of the expansion blocks* Only for the data presented in Tables 
XV and XVI and Figures 21 and 22 is the expansion block affixed to the 
downstream face of the swirler and measurements then taken at x/D = 0.0. 

4.1 Velocity Profiles From Radial Traverses 

Axial, radial and swirl velocity component data are tabulated in 
Tables III through VIII for radial traverses from the swirler centerline 
to the swirler exit radius. Data are presented for five values of swirl 
blade angle; zero (no swirler), zero (with swirler), 38, 45, 60, and 
70 deg. Corresponding velocity profile plots are shown in Figure 9 to 
14, with the profiles extending from the centerline to twice the exit 
radius (r/D =0.5 where D is the test section diameter used in associated 
studies). All velocities shown are normalized with respect to the 
swirler inlet uniform axial velocity, deduced independently from the 
pitot-static measurement upstream of the swirler. The outer ten data 
points are zero in each profile because the presence of the solid boun- 
dary of the swirler flange precluded measurements at these locations. 

The nonswirling case shown in Fig. 9 has a nearly-flat axial velo- 
city profile, as expected for the plain nozzle opening without the 
swirler installed. There is no measurable swirl velocity, and the ra- 
dial velocity is zero except for points very near the edge of the exit, 
where the flow begins to anticipate the abrupt expansion to twice the 
exit diameter. The second nonswirling case, see Figure 10, has the 
swirler installed with the blades set to cj) = 0 deg. The traverse was 
made midway between two blades and away from any of the hub supporting 
struts. Here again the axial profile is quite flat, with just a slight 
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inorease toward the hub. However, the velocity has Increased by nearly 
25 percent, because of the decrease 1n flow area with swirl er hub and 
vanes in place, In addition, the hub Induces a negative radial velocity 
across the entire annulus, overriding the tendency to anticipate the ex- 
pansion corner, The swirl velocity Is, as expected, negligible. 

The 38-degree blade-angle case In Figure 11 shows remnants of the 
flat Inlet profile over a small portion of the radius near the outside 
edge In both the axial and swirl profiles. The presence of the hub now 
constrains the three innermost points to zero, and the region between the 
hub and the flat portion In the axial and swirl profiles is approximately 
linear, The maximum axial velocity 1s 1,5 times the Inlet axial velocity 
because tlie flow area 1s decreased by the hub and also because centri = 
fugal effects have shifted the profile outward. The radial velocity has 
an irregular profile with a maximum value of one-half the inlet axial 
velocity. 

In the I = 45 degree case of Figure 12 the flat segments are no 
longer present and both axial and swirl profiles vary from zero at the 
hub to a maximum at or near the rim of the swirler in an almost linear 
fashion. The similar shape and magnitude of the profiles indicates that 
the turning angle is fairly uniform and only slightly less than 45 de- 
grees. The radial velocity is again irregular, but shows a step at r/D 
« 0.1 similar to that in the axial and swirl profilesj this 1s probably 
due to the central recirculation zone downstream beginning to slow down 
the flow upstream of it. 

Profiles ensuing from the case of ^ « 60 degrees, see Figure 13, 
all have a sharply peaked shape, with most of the flow leaving near the 
outer boundary. The radial component is considerably stronger, with a 



peak value nearly twice that of the reference velocity upstream of the 
swirler. The step in the 45 degree axial profile has now developed into 
reverse flow, indicating that the central recirculation zone now extends 
upstream past the exit plane* The reverse flow is accompanied by reduced 
swirl velocity and very low values of radial velocity. The positive 
axial velocity adjacent to the hub may be the result of a slight clear- 
ance between the blades and the hub, allowing air with greater axial 
momentum to pass through. 

Exit velocity profiles obtained for the strongest swirl case con- 
sidered (({) « 70 deg.) are shown in Figure 14. Almost all of the flow 
leaves the swirler at the outside edge. The maximum axial and swirl 
velocities are approximately 3 and 2.5 times the upstream reference 
values, respectively, and the velocity gradients across the profiles are 
quite large. The reverse flow in the center of the axial profile is 
stronger than in the 60-degree case and is now accompanied by negative 
or inward radial velocity. This suggests the possibility of a vortex 
ring structure occurring at the exit of the swirler under high-swirl 
conditions. The swirl velocity profile remains positive but shows a 
step corresponding to the outer boundary of the recirculation zone. 

4.2 Velocity Profiles from Azimuthal Traverses 

An indication of the azimuthal or 6-variation of axial, radial, 
and swirl velocities is now given for the same vane angle settings used 
in the radial traverses. The measurements were taken at a constant 
radial position of r/D = 0.179, which in most cases illustrates ade- 
quately the azimuthal flow variation. However, measurements at r/D = 
0.204 were necessary in the 45 ~ 70 degree case to get data more repre- 
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sentative of the main region of the flow. In addition, azimuthal tra- 
verse measurements were taken 0.109 D downstream (at x/D = 0.0, expansion 
corner with the 90-degree block installed) for - 70 degrees to inves- 
tigate further the upstream extent of the central recirculation zone. 
Radial profiles at this location for all degrees of swirl are already 
available (3). 

Measurements in each case span an angle of 48 degrees, somewhat more 
than the 36 degrees between successive blades. Data are tabulated in 
numerical form in Tables IX through XVI, and corresponding velocity pro- 
files are given in Figures 15 through 22. 

The variations in all normalized velocity components u, v, and w 
occur in approximately 36-degree cycles, coinciding with the blade 
spacing., The profiles all show significant variation With azimuthal 
position, except for those in or near recirculation zones where the W 
velocity component is dominant. These variations can be attributed to 
several causes, among them being blade stall from using flat blades at 
high angles of attack and wakes from blunt trailing edges. 

Figure 15 shows the azimuthal profile with the swirler installed, 
but with the vanes set to zero angle. The G - 0 degree position is di- 
rectly downstream of one of the swirl vanes, approximately 3 mm from 
the trailing edge at the r/d = 0.179 position. The velocity defect in 
the wake of the blade is clearly seen in the axial velocity profile, al- 
though the precise accuracy of these measurements is uncertain because 
of the velocity gradients across the width of the probe. The decreased 
u-velocity at the left side of the profile is caused by the presence of 
an upstream strut supporting the hub, located at 9 = +24 degrees. The 
radial velocity is uniformly negative indicating inflow over most of 
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the range, which agrees well with the results of the radial traverse 
shewn earlier in Figure 10. The radial velocity is positive only in the 
blade wake region. The swirl velocity, as expected, is effectively zero. 

Figure 16 presents the results of an azimuthal traverse for the ({) 

- 38 degrees Tow-swirl case. The measurement position at r/D = 0.179 is 
in the middle of the flat portion of the radial profile, as may be de- 
duced from observation of Figure 3. The 36-degree cyclic variation from 
one blade to the next is apparent in each of the profiles. The u and w 
profiles have a flat portion, apparently between blade wakes, with an 
average yaw angle of about 39 degrees. This confirms the assumption 
that the blade pitch/chord ratio of 0.68 is sufficient to adequately turn 
the flow. In fact, over the rest of the profile, the turning angle is 
even higher than the blade angle (j>. The radial velocity shows no flat 
region and varies the most of the three components. It is also quite 
large even at this low degree of swirl. 

In the case of ^ = 45 degrees, Figure 17 illustrates that the 36- 
degree cycle is not as clear, but nevertheless significant variation 
exists in all profiles. The radial component is nearly as large as the 
axial and swirl components in some places, and again exhibits the great- 
est variation with azimuthal position. 

For the 60-degree swirl case of Figure 18 variations with azimuthal 
position are again evident in all profiles. The variation is less than 
in the cases seen heretofore, r possibly because the main flow has shifted 
further outward under centrifugal effects and the measurement position 
is in a region of reduced velocity. 

This effect is even more notable in the c}) = 70 degrees profiles 
portrayed in Figure 19. The measurement position is now no longer in 
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the main exiting flow, but on the edge of the central recirculation 
zone. The axial velocity here Is effectively zero, although consider- 
able swirl and radial velocities are present. The radial velocity, it 
should be noted, is negative or inward towards the centerline. Azimu- 
thal variations are fairly small here, which is to be expected since 
the flow is mainly in the azimuthal direction. To get a more represen- 
tative sample of the exiting flow from the swirler with blades at 70 
degrees a traverse was made at the next outward radial station at r/D 
s 0,204. When the velocity profiles shown in Figure 20 are compared 
with those in the previous figure, the effects of extreiiie velocity gra- 
dients in the radial direction may be perceived. The accuracy of the 
radial velocity and pitch angle measurements may be suspect in the pre- 
sence of high radial velocity gradients, but the major features of the 
flow can still be assessed. In a radial distance of only 7.6 iiin, the 
axial velocity jumps from zero to over 12 m/s. In addition, the swirl 
velocity increases over 50 percent and the radial velocity changes sign. 
The 36-degree cyclic variation with blade spacing is again present in 
all profiles. 

To investigate further the complexities of the flow with swirl vane 
angle ^ * 70 degrees, azimuthal traverses were also made 3.25 cm down- 
stream of the location of measurements just discussed. Both radial 
locations, r/D ~ 0.179 and 0,204, were investigated at x/D - 0.0, This 
is the axial location of the expansion station in practice, (1,3,5, 

7-9) and the 90 degree expansion block was affixed to the downstream 
face of the swirler for these measurements. The profiles appear in 
Figures 21 and 22-, they may be compared with corresponding profiles 
from further upstream in Figures 19 and 20, respectively. It appears 
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from both sets of profiles that the recirculation zone has narrowed 
somewhat with the additional length before the expansion corner. At 
the inner radial position (r/0 » 0.179) of Figure 21, the axial velo- 
city is no longer zero. It is now positive, Indicating that the main 
exit flow has moved slightly further inward. The azimuthal variation 
is still quite small, however, suggesting that the damping influence 
of the recirculation zone is still in effect. At the outer radial posi- 
tion (r/0 = 0.204) of Figure 22 the axial and radial velocities are 
larger than at the upstream position, also implying that the outer high- 
velocity zone has moved further imyard. The azimuthal variation is 
again similar to that of the exit-plane position at the same radius. 

4.3 Calibration Sensitivity Verification 

Since minor variations occur from one probe calibration to the 
next, it was decided to check the sensitivity of the data reduction 
procedure to these variations. The case of swirl vane angle (}> = 70 
degrees was used, at x/D = -0,109 and r/D = 0.179. The most recent 
calibration provided the baseline values of the pitch and velocity 
coefficients, (5,7) which were then varied by increasing the magnitude 
of each value by ten percent. Three cases were tried: increased pitch 
coefficient with baseline velocity coefficient, increased velocity co- 
efficient with baseline pitch coefficient, and increased values of 
both coefficients. The percent difference in the output values of the 
velocity components is shown in Tables XVII through XIX for each of 
these three cases respectively. 

Referring to Table XVII, changing the pitch coefficient value is 
seen to affect the radial component the most, as expected. The change 
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in output stays below ten percent for all but three of the output values. 
For the case of increased velocity coefficient only, Table XVIM shows 
a quite uniform increase of less than five percent over all the values. 
This indicates a relatively predictable, low sensitivity response to 
changes in the calibration velocity coefficient. 

The final case, shown in Table XIX, indicates that increases in 
both coefficients tend to cancel each other for the radial velocity mea- 
surement, which was the most sensitive to pitch coefficient variation. 

The axial and swirl components increase somewhat, but all variations 
remain well below ten percent. This relative insensitivity to calibra- 
tion errors is satisfying but it should be noted if the coefficient 
changes are of opposite sign in the combined case, errors of greater 
than ten percent in the radial velocity measurements would probably 
ensue. 

4.4 Swirl Strength Comparison 

For comparison with the results of the idealized profile deriva- 
tions, swirl numbers S and S' were calculated from experimental data 
using Equations (1) and (6) with the turbulent stress terms omitted. 
Measured velocities and pressures from the radial traverses described 
in Section 4.1 were used, with appropriate numerical integration per- 
formed by the computer data reduction program described in Appendixes 

C and D. Since actual wall static pressure measurements were unavail- 

♦ 

able, the reference pressure was taken as the static pressure mea- 
surement at r/D = 0.230, the point nearest the outer edge of the 
swirler. The results are given in Table XX, showing the asymptotic 
behavior of the flat swirl vanes in producing strong swirl. Also shown 
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in Table XX is the ratio for each vane angle, taken from the 

mo mo 

measured radial traverse data. These ratios were used to compare the 
actual profiles with the idealized ones. 

Two comparisons were made to investigate the usefulness of the 
idealized profiles. In the first, swirl numbers from the measured pro- 
files were compared with those predicted by the Case I idealization. 

This was done by making the standard assumption that an ’’ideal'* flat- 
blade swirler (with an infinite number of infinitely thin blades) 
operating on a plug flow would produce flat exit profiles as shown in 
Figure 1, part (a). The flow turning angle would be everywhere equal 
to the vane angle and the ratio w^/u^ = F would be equal to tan <}>. 
Corresponding S and S' values for each vane angle are then found using 
Equations (13) and (15) or Figure 2 with F = tan <|>. The results for 
the four swirl vane angles used are shown in the left half of Table XXI. 
It is immediately apparent that the negative S values for c|) = 60 and 
70 degrees are based on values of F greater than the asymptotic value, 
and are physically unrealistic. The S values for ((> = 38 and 45 degrees 
are considerably higher than the measured values, while the S' values 
start close to the measured ones but diverge rapidly at high vane 
angles. This confirms the unsuitability of the Case I idealization 
for modeling flat-bladed swirler performance. 

The other comparison was done using the "most appropriate" idea- 
lized case, as judged by visual comparison of the profile shapes. The 
measured value of the ratio of maximum profile velocities from Table XX 
was used instead the tan cj) assumption, which has no theoretical basis 
for Cases II-V. Most appropriate cases were determined to be Case I for 
(j) = 38, Case III for <j> = 45, and Case V for (f) = 60 and 70 degrees. S 
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and S‘ values were then determined using Equations (13) and (15), (18) 
and (19), and (22) and (23). Results are shown in the right-hand side 
of Table XXI. Again we see considerable discrepancies between the actual 
and idealized values for both S and S’. Although use of Cases III and 
V gives a much better match for the higher swirl vane angles, the newer 
idealized profiles are still inappropriate for modeling actual swirler 
output. The disparities may be attributed to the presence of the central 
hub, the upstream extent of the central recirculation zone, and flat 
swirl -vane ineffectiveness at high angles of attack, with associated 
wakes and nonaxi symmetries. 



CHAPTER V 


CLOSURE 

5.1 Summary and Conclusions 

This study has investigated the perfoniiance characteristics of an 
axial vane-type swirler, used in combustor flowfield measurements and 
turbulence modeling research. A theoretical analysis of swirl numbers 
associated with several idealized exit velocity profiles is included, 
and values of the ratio of maximum swirl velocity to maximum axial ve- 
locity at different swirl numbers are tabulated for each case. Measure*^ 
ments of actual swirler exit velocity profiles were made for swirl vane 
angles ~ 0, 38, 45, 60, and 70 degrees using a five-hole pitot probe 
technique. The values of normalized velocity components are tabulated 
and plotted as part of the data base for the evaluation of flowfield 
prediction codes and turbulence models. 

Assumptions of flat axial and swirl profiles With radial velocity 
equal to zero were found to be progressively less realistic as the 
swirler blade angle increases. At low swirl strengths (ii 38), por- 
tions of the u and w profiles remain flat while the v-component is al- 
ready significant. At moderate swirl <{> - 45 degrees, approximately 
linear profiles of u and w with radius are found, with strong v velo- 
city. At stronger swirl <P = 60 degrees, even more Spiked profiles are 
seen with most of the flow leaving the swirler near its outer edge, 
and some reverse flow near the hub. At strong swirl 4-70 degrees. 
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the profiles are extreniely spiked with flow reversal , The central re*^ 
circulation zone extends upstream of the exit plane, almost to the 
swirler blades in high-swirl cases* Because of this recirculation and 
the presence of the hub, none of the idealizations considered could 
model actual swirl cases adequately. 

The flow-turning effectiveness of the flat blades was generally 
adequate for all vane angles tet.ted, However, the large variations of 
flow angles and velocities with radius made meaningful comparisons with 
two-dimensional cascade data impossible. Nonaxi symmetry was found in 
all swirl cases investigated. 

B.2 Recommendations for Further Work 

Other aspects of swirler performance not covered by this project 
include pressure drop across the swirler and the efficiency of swirl 
generation, It is recommended that these be investigated for the pre- 
sent swirler to allow comparison with values quoted by other swirl 
researchers . 

Development of idealized profiles accounting for annular flow 
and recirculation is another area in which further work is reconwended. 
This should include relating the ratios at maximum profile velocities 
to effective vane angles to allow prediction of swirler output for a 
given vane angle setting, 

Finally, it is suggested that an uncertainty analysis be done on 
the five-hole pitot technique to estimate the effects of turbulence 
intensity and velocity gradients on the accuracy of measurement re- 
sults. 

ORIGINAL PAGE 18 
OF POOR QUALITY 



ORIGINAL PAGE IS 
OF POOR QUALITY 


REFERENCES 


(1) Lilley, D. G., and Rhode, D. L., A Computer Code for swirling Tur- 

bulent Axi symmetric Recirculation Flows in Practical Isother- 
mal Combustor Geometries, NASA CR-3442, Feb. 1982. 

(2) Rhode, D. L., Lilley, D. G., and McLaughlin, D. K. , On the Predic- 

tion of Swirling Flowfields Found in Axi symmetric Combustor 
Geometries. ASME Journal of Fluids Engng., Vol, 104, 1982-, 
pp. 378-384. 

(3) Sander, G. F. , Annular Vane Swirler Performance. Proceedings, 

Thirteenth Southwestern Graduate Research Conference in 
Applied Mechanics, Norman, Oklahoma, April 16-17, 1982, pp. 
274-278,. 

\ 

(4) Abujelala, M. T., and Lilley, P. G., Confined Swirling Flow Predic- 

tions, Paper AIAA-83-0316, Reno, Nevada, Jan, 10-13, 1983. 

(5) Rhode, D. L., Lilley, D. G., and McLaughlin, D. K., Mean Flowfields 

in AxisymmetHc Combustor Geometries with Swirl, Paper AIAA 
82-0177, Orlando, Florida, Jan. 11-14, 1982. AIAA Journal, 

1983 (in press). 

(6) Lilley, D. G., Turbulent Combustor FI owfi eld Investigation. Paper 

in Combustion Fundamentals Research Conference, held at NASA 
Lewis Research Center, Cleveland, Ohio, Oct. 21-22, 1982, 
pp. 152-168. 

(7) Yoon, H. K., and Lilley, D, G., Five-Hole Pitot Probe Time-Mean 

Velocity Measurements in Confined Swirling Flows. Paper 
AIAA-83-0315, Reno, Nevada, January 1.0-13, 1983. 

(8) Janjua, S. I., McLaughlin, D. K., Jackson, T. W. , and Lilley, D. 6., 

Turbulence Measurements in a Confined Jet Using a Six Orienta- 
tion Hot-Wire Probe Technique, Paper AIAA 82-1262, Cleveland, 
Ohio, June 21-23, 1982. 

(9) Jackson, T. W. , and Lilley, D.G., Swirl Flow Turbulence Measure- 

ments Using a Single-Wire Technique. Paper AIAA-83-1202, 
Seattle, Wash., June 27-29, 1983. 

(10) Janjua, S. I., and McLaughlin, D. K., Turbulence Measurements in a 
Swirling Confined Jet Flowfield Using a Triple Hot-Wire Probe. 
Report DT-81 78-02, Dynamics Technology, Inc., Torrance, CA, 

Nov. 1982. 


30 


ORIGINAL PAGE S3 

OF POOR QUALITY 


31 


(11) Kerr, N. M., and Fraser, D., Swirl. Part I: Effect on Axi sym- 

metrical Turbulent Jets. Journal of the Inst, of Fuel, Vol . 
38, Dec. 1965, pp. 519-526, 

(12) Mathur, M. L. , and MacCallum, N. R. L., Swirling Air Jets Issuing 

from Vane Swirl ers. Part I; Free Jets; Part II: Enclosed 

Jets. Journal of the Inst, of Fuel, Vol, 40, May 1967, pp. 
214-245. 

(13) Chigier, N. A. and Chervinsky, A., Experimental Investigation of 

Swirling Vortex Motion in Jets. Journal of Applied Mechanics, 
Vol. 34, June 1967, pp. 443-451. 

(14) Beer, J. M., and Chigier, N. A., Combustion Aerodynamics. Applied 

Science Publishers, London, 1972. 

(15) Beltagui, S. A., and MacCallum, N, R. L., Aerodynamics of Vane- 

Swirled Flames in Furnaces. Journal of the Inst, of Fuel, 

Vol. 49, Dec. 1976, pp. 183-193. 

(16) Gupta, A. K., and Lilley, D. G. , Flowfield Modeling and Diagnos- 

tics. Abacus Press, Tunbridge Wells, England, 1983 (in 
press). 

(17) Beltagui, S. A., and MacCallum, N. R. L., The Modelling of Vane- 

Swirled Flames in Furnaces. Journal of the Inst, of Fuel, 

Vol. 40, Dec. 1976, pp. 193-200. 

(18) Morel, T., Comprehensive Design of Axisymmetric Wind Tunnel Con- 

tractions. Paper ASME 75-FE-17, Minneapolis, Minnesota, 

May 5-7, 1975. 

(19) Rhode, D. L., "Predictions and Measurements of Isothermal Flow- 

fields in Axisymmetric Combustor Geometries," Ph.D. Thesis, 
Dept, of Mechanical and Aerospace Engineering, Oklahoma 
State University, Dec. 1981. 

(20) Yoon, H. K., "Five-Hole Pitot Probe Time-Mean Velocity Measure- 

ments in Confined Swirling Flows," M.S. Thesis, Dept, of 
Mechanical and Aerospace Engineering, Oklahoma State 
University, July 1982. 

(21) Bryer, D. W., and Pankhurst, R. C., Pressure-Probe Methods for 

Determining Wind Speed and Flow Direction. Her Majesty's 
Stationery Office, London, 1971. 

(22) Lilley, D. G. , and Rhode, D. L. , "A Computer Code for Swirling 

Turbulent Axi syriwetric Recirculating Flows in Practical 
Isothermal Combustor Geometries," Report NASA CR-3442, 

NASA Lewis Research Center, Cleveland, Ohio, Feb. 1982. 






APPENDIX A 
TABLES 




original* page fS 

OF POOR QUALITV 


TABLE I 


RATIOS OF MAXIMUM SWIRL AND AXIAL VELOCITIES F-0 
OF IDEALIZED PROFILE CASES I - V, FOR COMMON 
VALUES OF SWIRL NUMBERS S AND S' 



om«Ai man rs 
OF poor? QUALITY 


TABLE I (Continued) 


$ 

0 

S' 

6 

0,10 

0,198 

0,10 

0.200 

0,g5 

0.47a 

0.25 

0,500 

0.50 

0.828 

0.50 

1.000 

0,75 

1.070 

0,75 

1.500 

LQQ 

2,236 

1.00 

2.000 

1.50 

1,442 

1,50 

3.000 

a. 00 

1.562 

a. 00 

4.000 


2,000 

\V 

v\T 


(b) Case II - Flat axial and 11 neat' sw1r1 profiles, 

W iLf />« 


C * w /u 
0 11)0 
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TABLE II 

SUMMARY OF OPERATING CONDITIONS 


(j> (degrees) 

FS (rpm) 

(m/s) 

Re^ X 10'^ 

0 

1950 

23.00 

2.22 

38 

2265 

13.30 

1.30 

45 

2600 

13.00 

1 .26 

60 

2800 

9.20 

0.90 

70 

2800 

5.52 

0.53 


* Abbreviations used are: 

^ Swirl vane angle 
FS Fan speed 

u. Spatial-mean swirler exit axial velocity, deduced from 
independent upstream measurement, excluding presence of 
the hub and swirler 

Re^ Swirler-exit Reynolds number based on and swirler 
diameter 
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TABLE XVII 

CALIBRATION SENSITIVITY COMPARISON, ACTUAL VS. 10% HIGHER 
PITCH COEFFICIENT ONLY 


Percent Difference 


K 

0 (deg.) 

u/Uin 

''/“in 

w/u,n 

1 

-24.0 

1.91 

-8.22 

1.91 

2 

-18.0 

0.80 

-10.23 

0. 80 

3 

-12.0 

0.27 

-11.43 

0.27 

4 

-6.0 

0.92 

-10.01 

0.92 

5 

0.0 

2.15 

-7.89 

2.15 

6 

6.0 

1.87 

-7.27 

2.87 

7 

12.0 

2.55 

-7.51 

2.55 

8 

18.0 

2.29 

-7.73 

2.29 

9 

24.0 

1.93 

-8.17 

1.93 
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TABLE XVIII 

CALIBRATION SENSITIVITY COMPARISON, ACTUAL VS. 10% HIGHER 
VELOCITY COEFFICIENT ONLY 


K 

0 (deg.) 

Percent Difference 
u/u,„ v/u,n w/u^„ 

1 

-24.0 

4.86 

4.86 

4.86 

2 

-18.0 

4.88 

4.88 

4.88 

3 

-12.0 

4.88 

4. 88 

4.88 

4 

-6.0 

4.88 

4.88 

4.88 

5 

0.0 

4.86 

4.86 

4.86 

6 

6.0 

4.88 

4.88 

4.88 

7 

12.0 

4.87 

4.87 

4.87 

8 

18.0 

4.87 

4.87 

4.87 

9 

24.0 

4.86 

4.86 

4.86 
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TABLE XIX 

CALIBRATION SENSITIVITY COMPARISON, ACTUAL VS. 10% HIGHER, 
BOTH PITCH AND VELOCITY COEFFICIENTS 


Percent Difference 


K 

6 (deg.) 

u/Uin 

''/“in 

w/uin 

1 

-24.0 

6.87 

-3.75 

6.87 

2 

o 

CO 

1 

5.72 

-5.85 

5.72 

3 

-12.0 

5,15 

-7.12 

5.15 

4 

-6.0 

5.84 

■ -5.62 

5.84 

5 

0.0 

7.12 

-3.41 

7.12 

6 

6.0 

7.88 

-2.75 

7.88 

7 

12.0 

7.54 

-3.01 

7.54 

8 

18.0 

7.27 

-3.25 

7.27 

9 

24.0 

6.90 

-3.70 

6.90 
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TABLE XX 

SWIRL NUMBERS S AND S' FROM RADIAL TRAVERSES 


4' 

S 

S' 

'^tno/%0 

38 

0.567 

0.559 

0.801 

45 

0.765 

0.718 

0.876 

60 

0.850 

0.759 

0.937 

70 

0.883 

0.750 

0.887 




TABLE XXI 





THEORETICAL SWIRL NUMBERS 
BY TWO METHODS 


<P 

Ideal 

S 

Case J 
S' 

Most Appropriate 
Case S 

Case 

S' 

38 

0.750 

0.521 

I 

0.786 

0.534 

45 

1.333 

0.667 

III 

1.137 

0.584 

60 

-2.309 

1.155 

V 

1.291 

0.625 

70 

-0.660 

1.832 

V 

1.066 

0.591 
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(a) Case I - Flat Axial and Swirl Profiles 



(b) Case II - Flat Axial and Linear Swirl Profiles 



(c) Case III - Linear Axial and Swirl Profiles 
Figure 1. Idealized Axial and Tangential Velocity Profile Cases 
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Figure 3. Photograph of Swirler - Upstream End 
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Figi' e 4, Photograph of Swirler - Downstream End 
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Figure 5. Diagram of Swirl er - Section and Downstream Viev/ 
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Figure 6. Swirl Vanes 
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Figure 11. Normalized Velocity Profiles From Radial 
Traverse, tj) = 38 deg. 
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Figure 12. Normalized Velocity Profiles From Radial 
Traverse, cj) = 45 deg. 
























Figure 13. Nonnalized Velocity Profiles 
From Azimuthal Traverse, (J> 
= 60 deg. at r/D = 0.179 
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Figure 22. Normalized Velocity Profiles 
From Azimuthal Traverse, (}) 
= 70 deg. at r/D ~ 0,204 
Measured 0.109 D Down- 
stream of Swirl er Exit 






APPENDIX C 

DESCRIPTION OF REVISIONS TO COMPUTER PROGRAM 
FOR FIVE-HOLE PITOT DATE REDUCTION 


APPENDIX C 


The data reduction program used for this project is a modification 
of a program written by Rhode (19) and described in some detail by Yoon 
(20). A brief overview of the entire program will be given, followed 
by a more detailed description of the major changes. 

1. Program Overview 

The reduction program consists of a main program, two function 
subprograms, and five subroutines. The main program first calls sub- 
routine INIT to initialize all array variables to zero, then reads in 
calibration data, control parameters, and the data to be reduced. The 
actual data reduction is done by repeated calls to the function SPLINE, 
which uses a cubic spline interpolation method to find pitch angle, 
velocity, and static pressure at each point from the calibration data. 
The function H and subroutines ABUILD and GAUSS are called from SPLINE 
as part of this process. 

Next a set of auxiliary calculations are performed. These include 
nondimensionalizing the output values, calculating momentum fluxes and 
swirl nuiflbers for radial traverses, and computing averages of the out- 
put quantities over successive one-blade cycles for azimuthal traverses. 
Finally, the primary output values are written into an unformatted 
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output data set for disk storage, and all output variables are printed 
out 1n standard format using the subroutines WRITE and PRINT. 

Changes were made to two sections of Rhodens original program; 
the main program and subroutine INIT. For brevity, only the changes 
to these sections will be considered In detail here. For information 
on the structure and function of the unmodified parts of the program, 
see Reference 20. 

2. Additions and Modifications 

The code’s new capabilities include calculation of static pressure 
at each location and reduction of either radial or azimuthal traverse 
data. For radial traverses, the code calculates axial flux of axial 
momentum (with and without statfc pressure contribution) and swirl 
numbers S and S’. For azimuthal traverses, it calculates averages of 
the output values u, v, w, and p - p^. In addition, substantial changes 
have been made In the way data is labeled, read in, and stored, in an 
effort to reduce storage requirements and make the code easier to use 
and understand. 

Static Pressure Calculation 

The static pressure is found using a method based on one described 
by Bryer and Pankhurst (21). The method uses the fact that the absolute 
pressure at any of the five holes in the probe tip can be expressed as 

p. = p . + K.q 

’^st V 

where p . is the local static pressure, K is an empirical coefficient 
which is a function of pitch angle 6 , q Is the local dynamic pressure 
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JspV , and the subscript 1 stands for any of the ports N, S, E, W, or C, 
Rearranging this and subtracting atmospheric: pressure from both sides, 
m obtain for the central pressure port 


Pst - Pabn = <Pc - Patm> ' 


(C.1) 


We now introduce the velocity coefficient, 

2 

which is already used in the code to determine total velocity magnitude. 
In accordance with standard practice, it is assumed that the velocity 
coefficients under calibration and measurement conditions are identical 
at a given pitch angle <S^, regardless of differences in fluid velocity. 


That is, = VC^^^meas 

Peal 


or 


’meas 


Tp 


c ■ Pw’«, ,cal ’ <Pc ■ Pw^«, ,meas 


This rearranges to 


’meas 


W^fi^imeas (C.2) 

Now, from Equation (C.l), taken at 5| under calibration conditions: 


"6^,cal 


Pc - Pst 


<Pc • Patni*«,.ca1 

_ I 1 


J(S^,cal 


xal 


since the static pressure equals atmospheric pressure in the free jet 
used for calibration. Substituting this and Equation (C.2) into 


Equation (C.l;. we get 
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’ ^atm^meas ” ” ^atm^<Si.ineas 


^*^C “ Patm^fi.j»cal 

• 

, -cal _ / V 

^cal 


(Pc - Pw>6,,ca1 


The calibration dynamic pressure cancels, and the remaining calibration 
pressures may be combined to form a dimensionless static pressure coef- 
ficient, 


SPC = 


(Pc - Patm> 

(Pc - Pw> 


which is determined as a function of <S from calibration data. This leads 
to the final expression for the gage static pressure at a location where 
the pitch angle is 6^: 


^•^st “ *^atm^6^ ,meas " " Patm^fi^.meas *■ " Pw^6^,meas 

This last expression is used directly in the code. The value of SPC is 
found by the same third-order spline interpolation technique used to 
find the pitch and velocity coefficients at each measurement location. 
(See lines 2690-2720 and line 3070 in the listing In Appendix D.) 


Radial and Azimuthal Capability 

The reduction of both radial and azimuthal traverses was implement 
ed by the addition of an integer flag in the input data to indicate 
which type of traverse is to be reduced. This flag, the variable 
KRADTR, is given a value of 1 for radial traverses and 0 for azimuthal 
Since this value is read in only once for the entire run, all 

c-g- 


oner 



86 


traverses to be reduced 1n a single run must be of the same type - 
either all radial or all azimuthal. 

Data for both traverse types is treated identically through 
Chapter I of the code, with the azimuth angles read in as values of 
radius, RINCHS. The major differences occur in Chapter II where the 
auxiliary calculations are performed. Depending on the value of 
KRADTR, radius values are nondimensionalized by the test section dia- 
meter or reset so that azimuth values remain in degrees. Next, the 
value of KRADTR is used to control branching to program segments which 
perform calculations unique to each traverse type, which are described 
in the next two sections. The last application of KRADTR is in Chapter 
HI, Output. Here again, it controls branching to ensure that only 
those output values appropriate to the traverse type being reduced are 
printed out. 

Radial Traverse Calculations 

When reducing data from radial traverses, the code automatically 
performs a simple numerical Integration procedure to find approximate 
values of mass flow rate and the momentum fluxes 6^, G , and G 
These values are then used to calculate the swirl numbers S and S' as 
defined in Chapter II. 

The integration procedure is effectively the same as that used by 
Rhode in his original reduction code, as well as in the STARPIC predic- 
tion code (22). However, the integration has been rewritten to cal- 
culate terms for the ring elements in a more straight-forward manner, 
and the central disk element has been added for completeness (lines 
3830 through 3880 of Appendix D). 


// 


m 

h\ tho absence of true static pressure taps In the rini of tlm 
swirleri the reference pressum p^ has been approximated by the measured 
static pressure at the nieasureuxsnt location nearest the wall of the 
swirler. This may Intitjduce an error, but the results will still be 
useful for comparing trends » 

Arlimithal Traverse Calculations 

For the azimuthal data, an averaging procedure Is used Instead of 
the integra'iilon routine. Since the data 1$ expected to be cyclic with 
a period of one blade width, averaging 1s perfo^ined over successive 
one«blade cycles. These successive averages may then be compared to 
chock deviation from cyclic behavior or averaged again to get a single 
representative value for each of the major output quantities. 

The code Is set up to handle traverses having six points over the 
width of one bladei for example, six-degme Increments for a ten«bladed 
swirl er. For other spacings the value of NREP (line 4470 of the code) 
must be changed. 

Since the reference pressure p^^^ for each vane angle setting 1s 
taken from a radial traverse at the exit plane, the value of p<* must be 
supplied by the user for azimuthal runs* This allows calculation of 
the pressure difference p - p^^ fjxjm azimuthal traverses for comparison 
with the values obtained from radial traverses. For those users not 
concerned with static pressure n\easuren«nts , the supplied reference 
pressure FREF may be omitted or set equal to atmospheric pressure. 

Mi seel 1 aneous Modif 1cat1 ons 

To make the code easier to use, all primary user inputs have been 
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se^^arated from the body of the code and incorporated into the block of 
input data, which is stored in a separate dataset. This minimizes the 
need to make changes in the body of the code, and reduces the memory 
space required to keep a record of all input data for each run. New 
headings were added to the input dataset to identify both the calibra- 
tion and measurement data, and additional variables are stored on disk 
for use by auxiliary programs which produce tables and profile plots. 
To improve readability of the code, all DO loops were indented and 
extensive comments were added. A listing of the reduction code with 
sample input and output appears in Appendix D. 


APPENDIX D 


LISTING OF FIVE-HOLE PITOT DATA REDUCTION 
PROGRAM WITH SAMPLE INPUT DATA 
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90 


ooouo 

00090 

oo;too 
00.110 
0011*0 
00130 
00140 
00 11^0 
oo:u-o 

00170 
001 BO 
00190 
00700 
001*10 
001*70 
007. -50 
OOP.40 
00011(0 
00260 
001*70 
001*00 
00 .'90 
00300 
00310 
00320 
00330 
00340 
003;ii0 
00360 
00370 
003B0 
00390 
00400 
004:1,0 
00420 
004:'50 
00440 
004 'liO 
00460 
00470 
OO4B0 
004‘'*0 

ooaoo 

003:1.0 
0031*0 
00330 
001J40 
001330 
00360 
001570 
00 3 BO 
00390 
00600 
00610 
00620 


Poor 


Page is* 
Quality 


D 

C ' 

C 

li 0 COMPIITKR P-ROGRAM F0F5 DATA REDUCTION OF FI.OE-HOLE PITOT 

C MEASUREMENTS IN TURDUUEMTr SWIRLlNOr RECIRCULATING l-XOW 

C IN COMDUSTOR OeOMETRlEB 

C 

C OERSION OF MARCH. 19B3 

C MOniFICATIDNB INCLUDE COMBINED RADIAL AND AZIMUTHAL. CAPA- 
C BILITY. REDUCTION OF STATIC PREGSULTE DATA. AND CALCULATION 
C OF MOMl-.NTUM FLUXES AND SWIRL NUMBERS FOR RADIAL 'PR'OPILES. 


0 

C BASED ON A PROGRAM BY- B. L. RHODE <PHD THESIS. OSU. 1981) 

C 

C 

0 G. F. SANDER 

C MECHANICAL AND AEROSPACE ENGINEERING 

C OKLAHOMA STATE UNIVERSITY 

C STILLWATER. OK 7407S 

C 
0 

si/ «i/ \f/ *1/ %t/ \l> W ^ ^ ^ iL> «i/ ij. ^ W W 

» t / ^ ' T * * * ^ ^ ^ ^ ^ ff • ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ■* ^ ^ ^ ^ 

C 

L‘. MAJOR FORTRAN VARIABLES IN MAIN PROGRAM (LISTED IN ORDER 

C OF FIRST OCCURRENCE IN THE PROGRAM)! 


IWRITE - LOGICAL FLAG FOR WRITING INTO OUTPUT DATASET (UNFORMATTED) 
BIAGMS - FLAG FOR DIAGNOSTIC OUTPUT 

,I:T - MAX no. OF TRAVERSES ALLOWED) DIMENSION VALUE IN SUBROUT INI:; S 

JT - MAX NO. OF POINTS ALLOWED PER TRAVERSE) ALSO DIMENSION VAl.UF 

HEDM ETC. - ALL VARIABLES STARTING WITH ■HED* ARE ALPHANUMERIC ARRAY'! 
FOR OUTPUT .HEADINGS 

MCAL “ NO. OF CALIBRATION DATA POINTS 

CI^ITCH - CALIBRATION PITCH COEFF. — (PN“PS)/(PC-PW) 

CDELIA - CAL. PITCH ANGLE STANDARD RANGE -5S TO •(•SB DEG* 

OVFI.CF “ CAL., VELOCITY COEFF. — (CAL, DYN. PRESS. )/(PC-PW) 

- CAL. STAl'IC PRESSURE COEFF. (PC"PA)/(PC-PW) 

.HEDID2 - CSFR HEADINGS I'D lUtNTlFY THE RIJN DEIHG REDUCED 

- INLET SIDEWALL EXPANBI.PN ANGLE 
“ SWIRL VANE ANCLE SETTING 

- INLET NOZZLE OR SWlRL.ER DIAMETER. DSMALL. IN ItNCHES 
" TEST SECTION DIAMETER. DLARGEr INCHES 

- INTEGER FLAG FOR TRAVERSE TYPE 1 FOR RADIAL. 0 FOR AZIM. 

- NO. OF TRAVERSES TO BE REniJCEB 

- MAX NO. OF POINTS IN ANY OF THE TRAVERSES PlIING REDUCED 

- AXIAL POSITION OF EACH TRAVERSE. INCITES 
■" NO. OF DATAPOINTS IN EACH TRAVERSE 


CT-STCF 

ITFDIDl 

ALPHA 

PHI 

ns INCH 

ril.lNCH 

KRADTR 

N.S I ATN 

MAXJPT 

XINCHS 

NDiTTA 

RDNPRS 

FREE 

T ANSPD 

TFLOW 

PATH 

B/OFF 

RINCHB 


INLET DYNAMIC PRESSURE (UPSTREAM OF 


REF. PRESS, 
FAN iipFEn. T 
TEMPERATURE 
ATMCSi'MFRIC 


USED TO CALC, PIUFF FOR 
PM 

OK A:l‘R IN TEST 
PR) SSURIT. TURR 


SWIRLER). TORR 
SWIRL. NUMBER. TORR 


SECTION. DEG. CELSIUS 


BETA /FRM-riFFSFr FOR V AW 
RADI AI. PCS. OF liATAPCnNT. 


ANGLE RFABINGS 
INCHES (THETA 


FOR AZIM. TRAVERSES) 



OU/itiO 
OOAi'.O 
00 A /l) 
OOf.00 
OOA«0 
00/00 
00/10 
oy/00 

00 /.‘JO 
00:^40 
00 >hi) 
00 /AO 

00 //o 

00 /BO 
OO/'/O 
00 BOO 

oortio 

OOB/O 

OOB.fO 

00U40 

ooniAi 

OOBAO 

QOB/0 

OOHBO 

oonyo 

00900 
00910 
009/0 
00930 
009.10 
ooyiio 
OOVAO 
009/0 
009BO 
00990 
01 000 
010 10 
oio:m 

01030 

01040 

01000 

OIOaO 

010/0 

0. 10150 
01090 
01100 
on 10 
01 1/0 
01130 
01140 
OllbO 
OllAO 

oim 
on BO 

01190 

oir’oo 

01/10 

01/:'0 

01 . no 

01. '40 
OU’bO 
01/ AO 

oi:‘/o 

01/00 
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1 ' ivl.l iO KAW MOllIt t(| "taw ANOM H la* 14 

I* KKNmb' » MIAB. OAtn 01 INOKIH mmin 1 Kl ‘.U, lUH ► lUKK 

a fOM'Biu ’ MtAU. vAt ta: tu mfmu = rwnn i iiibk 

0 BWU’A •• HIAB. VAlUl 01 I’t nm K iAlNtlBHUIO . lOKR 
C BeHOtl " INIU Nu.vn Ok OOlkHK RAOUlBf hlllKB 
0 kiAKiH - mn HimoN kAmuo* Hin ku 
fi X - AXIAt FtlBinoN 111 iKAVlkBbi MEUKB 
C t{ » KADtAl. t'tlHinON 01 UAl AnUNt» Ml It KB 

0 two - K AO 10 UBfc IMIRI l lUNi HK IH BPL INl INlf KPOl ATUIN KniU lNl 

i3 i'inu‘f . Rf uoftn piiT^f oom . lUk i-aph nATAtniiNT 

c lu.i lA “ Kitnjrni wmi ANon founb ba iNTikPiUATioN obinb 

G V'lnri • KKumn omnrm ootn. from iNtEkPOi atibn obinb oiita 

G porn - kroiin 11 uTATit* pkobb, oom . prom tniirpoi miun obino oh ia 

B RHO ^ lif-NBItY niK 1 AUH tRAVncBl * 1 RUM lOlAI -OAB I AM 

G HE1A - REDOPFD VALltl lOR PRtlBi YAW ANfU b r Obth 

1) VnVlAL ~ TUTAl VFinUlTY VEaiOR' NAUNTTUDEr N/B 

C O - AXIAl COMPONINI (IF OFlUCllYf N/B 

G V « RADiAR ttONP. OK Oil BOITY* M 'B 

{' W « TANtU-NTIAl (BWIRL) VH CICITY. M/S 

i: P - REDUCPH OAI III bp SlArtP PRt SSURla N/SU. M (OAOfcJ 

C XNH •* NUNWMENSrONAi AXIAL POSmUNf X/IH ARtlfe 

D BIN « INLET RH ERPNOf VPlOruV »CALC. IRBM RtlNPRSH M/B 

n MABMO * TNUT MASS FI BW RATb l ASSUMING UNIIORM AXIAl Vl pm’IIYn KB B 

G VTS1AR - NONWM, ruTAl OLl 00111 'lAUNtTOnia VTOTAI/UIN 

C USTAk “ NUNPIM» AXIAl ObtUFIlY> U ’UIN 

V VSTAR - NONOIMt RAWAl OPUJCms LVIITN 
V. USTAk “ NONniM. TANBPNTIAI VH , i W/OtN 

P F’STAR - NONWM. BlAn'* Pkf'SSURt r P RUNPRS 

C RNI.1 - NONDIM. RAWAl PUS.. R/IU ARBHJ ALSO IHEIA FDR A/IM. TRAOPRB! B 
0 nypii "> •Dfcl 'IA-Yf PUINT Stmr0‘ (I UK RAIUAL INTPORAIIONI I RBM BIARPU-i 
G IlYNP ” MibLTA- V» NORIH F’OINI* iSlM. TO HYPS) 

C PNS - 'UHAll NORTU BOLIlH* FROM S1ARFMF? USKI1 AS OH lA K FOR INTtOto. 

C PDIF F ~ PRbSB» Wn, F' * PRFF USPH TO CAI CUl ATP SWIRl NOMW R. N BIU H 

n ARPAl » AREA OF MBF Ht:HlNT AT Or NTFR OF JNTt CiRATlON RUmiN 

C FLOW - BUMHATION FDR HASS FlllW IHKUUGH RIND t lEMFNTB 

B WHOM » SUMMATION FOR ANGULAR HONFNTUM H UX 

U IIMUM « SUMMATtON FDR IiYNAHIt' AXIAl MOM. FU4IX iNPOl . PRPBH. it RM' 

V. tlMOMP - SUMMATION FOR AXIAl MUMINTUM FM1X> INiB . liTFSHURF Wfl , URM 

C AKEAJ =• AREA OP t AFH RING U f Ml NT> BO. M 

C MASS - INTFBRATFH MASS H DU RAIL. KD/S 

C IIMFAN ■* INIFORATEH MFAN AXIAL Old DPI TY. H/B 

V ANBMUM “ INTEGRA rpn AXIAL ROX OF AHGUl AR MOMENTUM. N- M 

C AXMOM - INI. AXIAL FLUX OF I1YNAMIC AXIAL MDH.» N INEGL. PK’IHS. llRi'n 

C AXMOMF' - INT. AXIAL FLUX OF AXIAL MONl NrUM^ N aNUL. P'Rl BviORE ITRH' 

C SPRIME “ SWIRl NIIMBFR OAU'. USING HYNAMin AXIAL MOMENTUM H UX 

F !i - SWIRL NlJMHbR CALF. 110100 FOIL AXlAl MOM. FLUX vINFl. I’REBB, v 

t UStAOO - AOERAOE of UBTAR OALUKB FOR AZIM. TRAV. » OVER ONt; liLAOt SPAl I 
e VSIAVG - AVG. OF VSIAR VALUL S 
0 USTAVG » AVG. OF" WGTAR VAl UPS 
C FI'FAVG •* AVG. OF PWFF VAl UPS 

U Visnos “ LAMINAR ASS. VlSCOSnY f Al f ULATl 0 FOR EACH IRAVl RSF f KO/Mlo. 

C KFWN *- INLET Kl Y'NaUiS WUMm iO OAl IS, USINO VISCUSITY FOR I ARM IRAVl 
C . 

0 

rUAPTLR 0 0 0 0 0 0 0 0 PRE I IMINARIES 0 0 Q 0 0 0 A 0 

-C. . . . . 

DIMENSION HEDH\*9> »HirmMNC9.1 .HF-DNMSEy UHF.OCMWCV) .HEDrMfYl9) • 
»HE0UC9>HlFnV<V).HFrai(91.HldlVrtV),HEDUBT(V1s 
+HF D VS T < 9 ) r HLDWS n 9 N . I ILDPS Ft 9 ) . HEOBE L 1 9 ) . HP DDE T 1 9 ) » 
iHL:DMMF{9^,j|t t’lNlvrvniU IiMlPi9) »HFDAM(9> r 

IHEDAX19) .HI DAXPCYY . HF DSPR(9> »HEDS(9) »*lb,DP(9) pHt.ilPDE 1 9) ,HEDKnu9) ♦ 
IllEDlBinS) .HLDUV.HIB) .HEBllBAt'?) jHt.DVSA(y) . HLDWBA<9) »IU JiPiiAlV) f 
IHLMIF'ANfy) iHFlVIFLaV) fHEDPAldV) »HEDRHtU9) iHEDVlS(9) ,Hl DCAl < 9) 

C 

COMHUN , 

I/CAL1D/CPI1UH(/A' rCPELTAlDA) .t:VEl CFCOA) MlPSTCF (HA) 


omoi^As: PAQl 

OF POOR QUAUry 


0_1L‘90 
01300 
013:10 
013120 
01330 
013^0 
013S0 
013A0 
01370 
01380 
01390 
Ol'lOO 
01410 
01420 
01430 
01440 
014 EiO 
014 AO 
01470 
01480 
01490 
OISOO 
OlSlO 
01520 
01530 
01S40 
01550 
015A0 
01570 
01580 
01590 
OlAOQ 
01610 
01620 
01630 
01640 
01650 
01660 
01670 
01680 
01690 
01700 
01710 
01720 
01730 
01740 
01750 
01760 
01770 
01780 
01790 
OIBOO 
018.10 
o;ui20 

01830 

01840 

01850 

01860 

01870 

01880 

01890 

01900 

01910 

01920 

01930 

01940 


C 

C-- 

C 

c 


/MEASUK/RUfclTA ( 8, 24 ) » Kf 'NHP8 ( 8r 24 ) » K’PCMPU < 8 r 24 ) » KPUI iPA 1 8 f 24 ) » 

NnATA ( t) > F MAX JPT f R1DNPR8 (lU , 

F0N8PI.1 < 8 ) » TFLOU ( 0 > i PATH ( 8 ) f hi f)FF ( H ) 

/GKQM/X ( fl ) F R < 24 ) f XNIi < H ) f KNU { 24 > f 1,'YPB < 24 ) » fJ YNP ( 24 > f 
BN 8(24) FNBTATNFXIHCHKaiT fRINCHB(24) 

/CALC/OTOTAL < 8 f24 ) f U ( 8 f 24 ) r V < 3 f 24 > f U ( 8 f 24 J f P ( 8 f 24 > r 

OTBTAR ( 8 F 24 ) F U8 T AK ( U f2.4 ^ f 9S TAR < 8 r 24 ) r WBTAR ( 8 * 24 ) » PBTAR i 8 f 24 > 
PTCHCF ( 8 F 24 ) f OKLCF C 8 f 2.4 ) f CELT A < 8 f 24 ) f BETA ( 8 f 24 ) . 

ANGMOH < 0 ) F UMEAN C B ) f HABB < 8 ) f MA8FL.C) ( 8 > f UlN < 8 ) f 
PE lIFF ( 8 F 24 > F PBTCF ( 8 f 24 > r AXM8M C 8 ) f AXMOMF' ( 8 > f 
SPRrME(8) f8(8) fREDIN(B) fPREF<8) FRHfJCa) f V.C8(::QB<8> f 
t IISTAVO ( 8 F 24 ) F OSTAye (8 F 24 ) f WSTAVG C 0 f 24 ) f PDF AOG < 8 f 24 ) 

<70UTPUT/BTC)RE<8) 

REAL MASSfMABFLO 
LOOrCAL IWRITEfMAGNS 

•BET TWR.lTt>'* TRUE. FOR WRITING SOLN. 8N DIBK BTORAGF? 

BET DIAGNB^.TRUE. TO ACTIVATE D1AGN0ST.IC WRITE STATPHENTS . 


IWRlTEa.TRUE. 

D1AGNS«. FALSE* 

ITaG 
JT=24 

•READ CHARACTER DATA FOR HEADINGS USED BY< BUDROUTINEH 
WRITE AND PRINT (ALSO CALISRATION HEADING) 

READ(5f205) HEDMfHEDUMNfHEDUfHEDVfHEDWf . 

HEDVT fHEDUST f 1-IEDVST f HEDWBTf HFDPST f HEODEI. f HEDDETf 
HEDNHS f HEDCMW f HEDCMA f HEDMMF r HEDMI V » HEDM.i;P f HEDAM i 
HEDAXf HEDAXP » HEDSPR f HEUS f HEDP f HEDPDF f HEDRED f 
HEDFAN F HEDTFL f HEDPAT f HFDRHO f HEDVIB f 
HEDLISA f HE'DVBA f HEDWBA f hi: DPD a f I IFDCAL 
205 FORMAT (9A4) 


—INITIALIZE VAR;i.ABLEa TO ZERO 

CALL INIT 

“——READ FiVE“HOLE PITOT CALIBRATION DATA 
NCAL«25 

* DO 10 I^^IfNCAL 

READ < 5 F 210 ) CPITCH < I > f CDELTA ( I If CVELCP < I ) f CPSTCF < I ) 

10 CONTINUE 
210 F0RHAt<4F10.S) 

IF(DIAGNS) WRirE(6F400) <CP1.TCH( I ) f I>=a f25) 

IF(DIAGNB) WR1TE<6f 400) (CtlELTAl 1) fI“-1f2S) 

IF<DIAGNB> WR,ITE(6f 400) iDVELCF (X) f-I^^^IfPS) 

IF(DIAGNS) WRITE(6f400) (CPBTCF(I) fT«1 f25) 

400 F0RMAT(///rlXFl3(F8,4FlX)F7/F5XFl2<F8.4) ) 

—READ USER HEADINGOf GEOMETRIC AND, CONTROL PARAMETERS APPLYING 
TO ENTIRE REDUCTION RUN 

REALU5f215) HEDID1fHEDID2 

215 F(.)RMATC1BA4) 

RE.AD(5*216) ALPHAj PHI f DBINCHf DLINCH 

216 FCJRHAT<4F10.5) 

READ (5 F 21 7) KRADTRfNSTATNfMAXJPT 

217 F0RMAT(3I10) 

•READ EXPERIMENT PARAMETTRB SPECIFIC TO EACH TRAVERSE f THEN 
ACTUAL MEASUREMENT DATA TN TRAVERSE 
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ORiatNAi: PA«3E 

OF POOR QUALITY 


lUViiO Itt) .^0 I- i»NH(ArN 

01 RfcADCfi. WO) XlNl H'iU ? »NllAJAC 1) f KllNfUim: ) rfW I ( i i 

OlV/0 Rl-ACK‘i»?Ui) ^A^OJ^•|f(J )»TFLUW<J)H••Ar^{<, U^I^/nf f iT> 

oiv!?o .ipin-RriATAn) 

Ol'^VO DO 20 J -irJt')U 

O'HHHi Rt AD C 0 y 220 ) lU IHM !ti < J ) y Rft( 7 A (I y . I > . Kl -MHl H U » Hi KMJHPW vX y J ) y 

020 iO #■ RI'OMPAayJi 

02020 20 UONTINUR 

020AO M) liUNTINlIh: 

02ai0 C 

l>2»r.O C»‘ - — . CPNyFRT X'S AND K'S FROM INCHtsI JD MtIbRO 

02060 t: 

o;! 0/ 0 RBMAU -nBINnH'KO , 02'-, 4/2 , 0 

02000 • ro.ARfit -DLINIJH’HO *0224/2 yO 

<)20*/0 no as I ■XyNSTATN 

02 J 00 X ( I ) -XJ MHI 10 < I ) JftO » 02S4 

02110 JPTS'NDATArr) 

02120 DC) 32 jH^yJPTB 

Of’jaO R<J)->RXNCHO( J)*A*02S4 

02140 32 nONTINUE 

021 SO 315 CONTXNtJP 

021A0 220 l■^)lv'MAT(!»^• 10*13) 

Oi>l/() 230 f DRMATaf-lOtS'y m0y2H0.S) . 

02,100 IJ'^rUAONP) WRITl-(Oy4 '0i <NDA1 Ai I. > . X- 1. » NB CA IN) 

02190 XPCDIAGNB) WRXTE(A»4D0 ) <Xa ) yX*'J yNSTATM) 

02200 IfCrrJAOW!) WRITP( /.ySt'O? CR(.)) » J*.l >.Jpt<:i) • 

02210 DU S>7 .t^"l»NHTArM 

02220 XF'(DIAC!NS>) URITr- icjySPO) th'Dl-TAC !».))» I-I, y.lPTU) 

02, ’■•’30 TF < II :i: ABNS > Wk t T 1- ( 4 y hOO ) ( RPNHPB J •• !, » .IP FS ) 

02240 I'PCDlAONfi) URITt- CAySOO) (PPf)MPlP T y „)) y ..)• I. y Jl * 1 U ) 

02'’S0 IF C m ABNB i WRIT! ( 6 » SOO ) ( RfCMPA •!».))» J- 1 * JP I S > 

O'APAO ,r/ ruNfifJur.:: 

02.2/0 4S0 f-0RMAn/y40My:UFe*4flX)) 

O22U0 420 f*0RMAT(///y40XylCH)y.lX) ) 

02,290 500 FaRHAT(///y20Xyi0<FB,4) > 

02300 V. 

02310 t'HAPfPR 1 1 :|, 1 1 data REDUl I ).C)N t 1. 1. 1 1 1 

02.520 n 

02330 C--» -»~‘“~”»CALq PXCHCF AND .TN rERPni ATE F'-R DPI, TA PROM 

02340 € PITOT CALIBRATION CUROE 

023S0 C 

02340 rrUD' 0 

023/0 DO SO I";LyNHiTATN 

023(10 ‘ JPTO--‘NDATAa) 

02390 DP 40 J--lyJPTS 

02.400 XFC(RPCMPWay.i)*F0,0«0)*AND.(RPNMPSn:y.l).EO*0.0)) 00 TO 30 

024 10 PICHCF 1 1 y .,) > ^^RPNMPS (t y J ) / ( RPCNPW <1,J > . E-6 ) 

02420 IPC (PXCHCI- (I* .0 ,l!t.2.S44)..DR* (PICHCF (I» J) ,LT*-3*7A9) ) 00 TO 30 

02430 IF(TDID *FC1* Oi DPI TAay.D^-SPLlNpCCP.tTCHy 

02440 II t7i!7„TAyHi:AI yPXUllCF 

024S0 IF ( (DID *01* 0) D1 LT AU y J ) -BP ( CPH CH y CDEt.TA y 

02440 II NUAI .y| ’.rPHCFn'y .J) ) 

024/0 J.DlD-7 

02400 6(1 10 40 

024V0 .50 CONTINIJF 

025(H) DFLTAay J) ^0*0 

Ol'tno WRX (7.(4^050) '( ..I 

02520 OSO F(IRNAT(20Xy'PiLTICF 16 OUT OF RANGE OF CAL t BRA'/ 1 UN AT 1 - 

02530 I 'yX3y' AND 

02S40 40 CaNtXNUE 

02SS0 50 CONT I' Ml if; 

02S40 C 

02S/0 L'--—-—— INTERPOLATE KIR MELCF AND PSICF FROM PITOT CALIBRATION DATA 

O'."'i00 0 

02S91' nUD 0 

02/00 DO 00 I^-lyMOTATN 
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OI'AIO 

02'61J0 

0:’>(f)30 

OS6M 

02660 

026VO 

02600 

026V0 

02700 

027.10 
02/20 
OV/IH) 
02760 
02/00 
02/60 
02. /‘/O 
02/B0 
Od/90 
02800 
O'20'JLO 
02820 
020a0 
028-1^0 
026ti0 
020A0 
02870 
02880 
02890 
02900 
029:10 
02920 
02930 
02940 
02950 
0297)0 
02970 
02980 
02990 
03000 

030.10 
03020 
03030 
03040 
03050 
0307)0 
03070 
03080 
03090 
03100 
03110 
03;L20 
03130 
03:140 
03150 
03:!, 7.0 
03:1.70 
03.1 80 
03190 
03200 
03210 
03220 
03230 
03240 
03250 
0327.0 


ommm ims o 

OF POOR QUALITY 


.jfTii NHOTOU) ■ 

no 70 j.ifJFro 

IF ( ( lO-DHPW ( I » J ) , EO , 0 . 0 ) . OMn . ( KF NMfr. ^ 1 , .) ) . HI . n , 0 \ ) ou i .) ^ . 

iF((ons<riruo(if3m .or. nii.o) no m 

IF<nHD 0> Vtl.CFa#J> OPLiNt U:r«l l. I7i* 

J OOFLCF^NCOl, iHEUTAClf J> > 

XF(IEiIi;i .GT, Oi OElCIFa rji. 8MCDFLIA»f;yf.l Of » 

I NnOLfHELTAarJ)) 

IF<ini£i .50. 0) H;Tt;F(T»J)-SfLXNF(Crif 1-TA* 

II CP8TL‘FrNC0l »ElbLTAa> JH 

xranxii .or, o) psrcp(x»j)i:.8F(ont,uft.i;pf5ri:Fr 

I NUAl»nFl1A(Ir..n) 

. i PIP-1 

OU TO 70 
AH nOHTlNUE 

05LCF CI, .1) -0,0 
pBiapaf„i)"0.o 
WkXTF<A»090) InJ 

890 F0RhAT<20X» 'PU, lA XB nUT OF PAN15E OF nAUPF'A I i (M In h, 
f At Xc'i'’H3>' AND J»'fX3) 


70 CON’nNUF 
HO i;,0HTiNur:: 

8 

no 85 X==^l.N8TnTH 

iFdiTAKMS) ukirrf^.'iOO) <PinHCFa».n-.!-i .jf tb) 

I F < ii 1 AHN8 ) won h ( 6 . ",)00 ) i PFFTA 1 1 , J j » V- 1 f ..IP I .? ) 
XF<B1A8HB) WRITF<or500) < OEI CF <. X » J H J ■ 1 » Jl' 

IF ( it] A8MB ) OKI 1 1; ( 7, » 500 ) ( PBTX'.F X .X » J H . vIP 1 8 ) 

85 CONTINUE 
U 

-CAU; HAGNITUliF 8!- TOTAL MEAN VFL 08ITY VtCn'OR!. 

0- U* Oy R W COHPDME^NtB, AND BTAt-U: Pin"88llRl 
C 


PI-3. 14159 
no 100 X>^lyNBTATH 

RHO ( t ) wPATH ( :i: ) 1 ( '1,33 . 33 ) / >, ->86 , 9 At (. TFU'iW ' X ) 42 51 .15)) 

jpTa-NHArAa) 

i;iO 90 0-l»JPTB 

i(KTA(.T» J)«37)0.4'&ZUFFC.l ) •RBE'T A ( ,H J > 

.l'F< (RPOMPWa» J) .EQ.0.0) iANii. ' RPNHPH- U .1) i- JUO.O ' Mtp I t * o.O 
OTDTAI (.t r .) ) - HURT ( APS C2 . 0/RHU 1 1 )»UFUF FRPPM W ■ 1 » • i • M..*T , V . ) 

I.KX. j>55vraTAL.af J) t i2m(ni!,i ■ia(i.j’'vim 7:i»),,o> k 
. # C0S(BL“TA<.X..l).TPI7:l,8(’ 0) 

Va*4>»VT0TAI •* S I N ( OFU A (I » J ) » P I / U.0 » 0 ) 

Ma.j).-yTC)TALUy.i) K (..‘Ub‘n:ife;i ia( i ...u^pi/ipo.O) * 

# SXNUltn AX j: ) J>*PH l80.O) 

P ( X y .,) ) » 7 ICPCMI -A 1 i: » J ) "4-'8 If r,F' 7 1 y J ) *Pi •CMF M ( I j J > ) * 1.13 . 33 
90 PONTlHUf- 
100 I2.INTXNU1:: 

If- (DTAONS) WRnFXA-nOO) (omiAt. 7 I ,,.J) D.i -J . jpth:- 

IFtnXAONH) lJRXTL7A,-500)7LhH.))y l 1» l!-Tvn 
iFaHARNB) WRM:T!.:i 7,»!i07)) 70(:i:,.)H.i 
IF (PXAtiNB) WF;I Tr 7'>.50O, (Uic:r» J) . ( I s IpJ B - 
|;FU;iJ:A 8N!;!) WRI t I- < Ay50o;» f pf ) * n - ,1 ) i , IPIB > 

C 

CHAPri R 2 2 2 2 2 2 7UI*<Il. f AFX'f F.iXLCULAl lOMS 2 2 2 2 2 

0 

C— --NONlUMeNBIUMALIii- I.FMBIHB AND 7'H I'lCX Hl-B 

c: 

DD 150 ;t-^l»NSTATN 

XHnXl)' 7( 1. )/(2,01 |.'L.aR8F) 

OPTB-NHATAa) 

UIN f n«<HUR r ( 2 . 0 7|:-HH < \ ) .-^siyliMPF ;» a>* 1 3.1 . 3.-S ) } 4 7 7, ,312 ■ 5 , 938) U2 
MABFt.()7.T) P14h’Fl'l(l)ttirN(n4:R5fF)l Ut2 
tiC) 140 JHy.jPTii 

OTBI'ARXXyJ) OUlliAL 7.i:».J>/llXN< I ) 
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oRismAL page JS 
OF POOR QUALITY 


0 .?LV 0 

i>Mmo 

032V0 

03300 

03310 

033JH) 

03330 

033-10 

033b0 

03360 

03370 


UHrARarJ) ) 

MHiAi{a»j’ ■ya.jj/uiNU) 

WSTAR ( I f ,1 J W { f r I ) /IJIN 1 1 ^ 

PHI OR (I f 3 ) P ( 1 f 3 ) A RONf iiH (1 ) * 1 33 , 33 ) 

MO euNriMun: 

11)0 tniHnNUF 

Jf-aHA'iNr}) URirt:(/)M;'.0) (UlNdlM .IfN'II'ATN) 

I F' ( II I AliNB > wra T t“ ( 6 r 4B0 ) i M ABFM) ( i ) f t- 1’ » NtJI A I N ) 

no 160 3^;MMAX3Pr 

RNIi ( 3) s<K < 3 > / (5> , 0)f'B't,AF';0E ) 

IF f isRADTR . EQ . 0 ) RNO ( 3 ) r^RlNCHB < 3 ) 


03300 lF(KKAnTR,EO,0) R(3) "RINi:iiO(3) 

03*370 160 aUNTlDUE 

03400 t; * 

03410 IF(KRAriTR)Ea.O) (it) TO 1'3S 

034P0 tl 

03430 f> “FOR RAniAl. PROFILES} NUMBRIOAI.. INTBBRATION TO CALO. r.A ?U 


03440 C 

034!i<> t: 

03460 a 

0*3470 0 

03.4 BO C 

03470 0 

03B00 

03E10 

03S70 

03530 

03540 

03550 C 

03U60 

03570 

03580 

03590 

03600 

03610 C 

0367,0 

03630 

03640 

03650 

03660 

03670 

03680 

036VO 

03*700 8 

03710 C 

03770 D 

03730 G 

03740 C 

03750 

03760 

03770 

03780 

')3790 

03 BOO (. 

038:1 0 

03870 

03830 

03840 

03850 

03860 

03870 

03880 

03B90 

03900 

03910 

03970 


FLGW AND MOMENTUM FLUXES FOR 8UIRL NUMWii 


FOR PROFILES AT AND UPSTREAM OF EXPANSION COFiNFR*- RSft/iM 
IS USED IN EXPRFSSIUNS TOR DYNP AND UMEAN? HOWNBTRLAM Ui 
f-XF»ANSIONx RLAROE IS USED, 

» 

no 130 .t'-^l^NSTATN 
JPTS.’aN0ATAa> 

JPTSM1«JPTS"1 

PREF(I>.-P<Ir-JPTS5 

riYPS<15~0.0 

IF < XINGHS <1 ) . OT . 0 . 0 ) SO W 107 
DYNP < JPTS ) «■<! ♦ 0)T ( RSHALL-R ( JPTS ) ) 

GO TO 108 

107 nYNP(JPTS) -2*0*<RI, AR6F"-R(JPTS> ) 

108 GONTINUE 

no no ..)-^i.»upTSMi 

nYNP(JOi)xRtJ'M)H;(..J)' 

nYPS(„H'l>--nYNP(3) 

110 CONTINUE 

no US 0“if..iPTS 

SNS ( J ) ")0 ) &'!({ < 1.1 YNf ( J ) 1 IiYPS ( J ) ) 

PWFFdfJI-Par.D-PKEFCI) 

115 GONTINUE* 

-“-INNER 3 (HUB) VALUES OF FWI-F AlcE SET lU 2!.l?0 !• Ol-; 5W;n:i EP 
EXIT-PLANE PR0IMLES5 POFil DOWNSTRE AM PROf'Ii„l- 8 » AGTUAi, 'AUJU) 
ARE U?!EIi 

iFCXINGMSa ? ,8S' (iO TO 116 

PDIFF<Ir:t) 

pnil-r- vUE) 0, 

P01FF*(U3) "0* 

116 GONTINUE 

IFaitAON',) WRIIECUSOO) • DVNF ( j ) r-I MJI'TS-) 

IF’(ni'AUN5‘) I'lRl i F' • , • > . ' 'ifs > < ( J ) if J I • i| ' I ' i 5 

ARFA1- Id tSHE< DYE * 

ARSUM^AIILAI 

FL()W' IOin(X).TU< :i * I .DAFiFAl 
WMOH- U ( X y 1 ) YR( 1 > ) /4 » IFUIW 
UMOfi U(Tfl)tn,riW 

UMOMP (Fli(T£I) MJ(J'v D*Y7.iPni:ro- {I vl))YAf‘:EA( 

IFUiIAON?’) t)II-:i FECAD’OTO ) APEAi -)'iRH1JMjI l,OW>f(' 
no :l'’0 JD'r.lF'TS 

AFOiAJ-))’',YPlYF'< J)^:FiN*H,n 
ARGUN- M'WiUMFAREAJ 


ORIGINAI mGi IS 

OF POOR QUALITY 
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0.JV30 

Qd'Hky 

Qi9h6 

0AV6Q 

OJV'/'O 

0^000 

Q4UjO 

040.*a 

a4<»i5» 

04040 

04<»h0 

04O/.0 

040 <^0 
04000 
0405*0 
04X00 
04U0 
<)4iyo 
04130 
04140 
041Ii0 

041 AO 
04170 
04100 
04190 
04300 
04310 
0433O 
04730 
04240 
04a?i0 
042A0 
041*70 
042f)0 
04390 
04300 
04310 
04330 
04330 
04340 
043ij0 
043AO 
04 J70 
04 300 
04390 
04400 
04410 
04 470 
04430 
04440 
0441(0 
044A0 
044/0 
04400 
04490 
(54001) 
04010 
041520 
04!(30 
04’’(10 
04!4(O 
04:‘AO 
04 '770 
040H0 


FI nw t-UDWlRHlh Jl*uaiJl *F)KI rt.i 
UHllM'Umirt H<HU 1 1 ) liU ( 1 1 .1 » litof AKl rt.l 
UHOHP (HHlH 1 "KUU » f f 1 , .1)) *rtRt 0 i 

UMOH ■ RHiJi 1 ) n’j( 1 1 J) 'fW ^ 1 » .1 >13v*f .1 » 

IFCDIAtmiil Wktn,< A 1-70401 M(tn3fOKUUHfrUJl''"WH(Ut»f)i)»Sfi 
130 CONTINUE 

■ MrtSlUD-ROW 

a 


JFlXJNfiHSU ) ,15T.0.01 00 I I I.*;* 

UHtYiNCI) Mftoon «/^h;HO{n*. 
on TO 133 

u*:* uMroNarrtooaa )/<mhh i RriiKm ARuriR(3 > 

123 aONTJNUr 
0 ■ 

ONOHOMm-sWNOH 

AXMOMCn^UilOM 

AXMOHpnjwummp 

IFfOrOONS) WRITF^AdSOGO) IIMFAN(I) »MrtBinrjfONliMRn 1*. 
I OXMOHP(l) 

f: 


3030 rnRMATC/74Xf 'ARFA.JRSX* "ARiiUM' »f(>:» 'FLOW * »*Xr ' WHOM , . - 
i ' UMOM ' » AX » ' UMOMI ■ V/ ' ' r A4l 0 , 3 1 

3040 rOi^HATR •♦aL10.3> 

1*050 fflRMA1(/l4Xi-'OMrAW-»l,Xr*HAOaRAX» -ANOMOH'MXy AXNOh'' 

A 5XRAXM0M!-"//m,5L10.3J 

0 



130 



OPRILiP ( 1 ) -^-ANOMUM (I ) ^ i AXMOM f 1 ) HR31HAI I ' 
OCX > --AHOHOH 1 1 1 / ( nXhOdP ? I ) IfROMOI k ) 


CON riNue 

1 f { n.i AONB 1 wRi ri- ( A Mino ) 
If (DIAONO) WK'ITfc (6,450) 
I r- ( D 1 ABNB ) WR I Tk t h » 450 ' 
JF< OIAONS) WRITF: (A » 400 ? 
IF ( n .TARNS ) WRI Tl- ( 6 , 4G0 ) 
TF'COTAONS) WR.nt‘ (/>,40(') 
IT ( I)I AGNS ) WRTTf ( A , 450 ) 
CONTINUE 


(UMLAN(T)»I".l,NSrATN3 
(MASSCI) »I---l»NSrAfN) 
(ANOMPMU',! RNOrAfM) 
(AXMOMvl R l' -l»NS rA rN) 
(AXHOMIHI RT l,Nt!TArN> 
(SPRIMEd),! l,N'SrATN) 
(SCDrl-'lrNSTArH* 


IF (KRAOTR(EO.l) RO TO IHO 


V . — FOR A7IMUTHAI. TRAOERSr-Sl CAI, C» POl fF ■'^(P-PRl'F) 
D VAI.ULT OF PREF(I). 

c: 

. DO 17a .X-J,NSTATN 
UPTatiNDATAa) 

DO 177 Jc-UrJPTS 

PDTFF( I » J ) "p ( 1 , ,J ) »PREF ( I ) *133 . 33 


U'lTNO »;rj| ' 


177 CONTINUE 
170 criNTlMUF 


lL“— -CALC, AUPRAOirUAMIFC FOR A/IMimiAl rRAOKRSfC NRLl- iti ni 

P PO.tNTS IN RFfTAf INIL CfnU;. AtCROiFi ONI” Dl AlJti « NA‘*L IS f)lL -ii 
P AUFRASES POSSll'UI nONrAlNINtT NRIP '(WirCUriWE POINIVL 

r\ 

NPfP'-A 

DU IBO r IfNCfAm 
NAUrL NDATAn ) ■ NCfPfl 
DO 175 K~t,NiTyf. 

NAUfND.KfNRrP I 
USON 0* 

VSUM 0. 

WSUH-0. 

PSOM -0 . 

DC) 174 ,i KrNA‘»fND 
USUN“Ui!UMI!f iTARClfJ) 

OSUM VSUMfM!rARa».l) 
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OAU'fO 

0^600 

O-iAlO 

04a;‘0 

0*1 A 

04A4U 

O-IASift 

^4/.A0 

04A/0 

04/4f«0 

OAiiVO 

04;0v) 

04 7 1 a 

047;»0 

04VAO 

04740 

04/fiO 

04/40 

04 7/0 

04/UO 

04770 

04B00 

04U10 

04}1?0 

04iJ30 

04H40 

0401jO 

04B40 

04(370 

04Hfi0 

O4itV0 

04900 

049t0 

0491*0 

*'^>.'30 

04 700 
04 9 AO 
04970 
04980 
04990 
OSOOO 
onolo 
orjooo 

0(7070 

o:j')4o 

•'(jOAo 
OfiO •’0 
oi.oiio 
0(| 09o 

034 100 

001 to 

OSK’O 
Oil I io 
00140 

o:>i!.o 
Ob I AO 
Ob, 170 
051 BO 
OJll 90 
o:;..'0(i 
0 

0*j7=«i 

0‘t.'4«) 


oRiGif^Ai: mQK m 
OF POOR QUALITY 


WtiUM WHUMHilHlAf'a U.li 
f'.um 1 arj) 

1/4 tumriNiJi 

U£tTAVSi( J»M umm -fJi.l P 
M£iTAVHU»Kl OHUM» faiH' 

WOlAVfiUrK) WSU3UHK1F 
FrUAVUvlffs) PHim-NIO'T 
1/5 coNrimiE 
13, ‘0 nnNtiNiif 
i; 

C -BAlCULATt yisoosm ANA INIE! KE'fNJlf US NUNW-R (lUlTH li ./‘i f 1 t . 

e 

0“» -oisfinsnY rirmula fr'Om i.an r j-oskah* airfiani: ai Miimhuib •. 

8 * 55 FFRt tIFiMANCFf P.42* 

U 

00 1A2. I«1»NS(ATN 

DEN8H rFU1WaH2/A,iriH10*4 
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The following listing is of a dataset containing the input data 
for the reduction code. The two datasets are submitted together as a 
single batch job; they are merged by the computer before execution. 
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Output generated by the reduction code using the example data given 


above appears on the following pages 
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